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Abstract 

We present a new finite-size scaling method for the random walks (RW) 
superseeding a previously widely used renormalization group approach, which 
is shown here to be inconsistent. The method is valid in any dimension and 
is based on the exact solution for the two-point correlation function and on 
finite size scaling. As an application, the phase diagram is derived for random 
walks with a surface-bulk interaction where the system has either a surface 
or a defect. Possible extensions to disordered systems are also discussed. 
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I. INTRODUCTION 



Renormalization Group (RG) has been a cornerstone in the analytical evaluation of the 
critical exponents of various statistical models in the past two decades [|l],|2]|| . Its application 
to self-avoiding walks (SAW) was a natural consequence of the formulation of SAW as n — > 
limit of the n- vector model whose two-point correlation function yields, in this limit, the 
generating function of the random walk with the excluded volume effect [II] . 

An extension of this approach to the case without self-avoidance was carried out by 
Family and Gould (FG) JjJ in the absence of disorder, and by Sahimi and Jerauld (SJ) |7J 
and by Gould and Kohin (GK) |§ in the disordered case, where the disorder was mimicked 
by a site (or bond) percolation cluster. The results in the latter cases were found to be in very 
good agreement with those of numerical simulations §J of random walks on a percolating 
cluster on the square lattice. However unlike GK who used the kinetic rule for the random 
walk (more precisely they solved the model which is nicknamed myopic ant) fL0| , SJ used 
the static recipe which is now recognized to correspond to the so-called ideal chain model 
TT| r p~2| . p~3| , pr^JT5|1 . Recent analysis p^^2|J13| , p^JT5| with various different approaches have 



shown that the ideal chain does not belong to the same universality class as the kinetic 



walks; rather it is equivalent to the random walk in a trapping environment [H],P|| • 111 view 
of this, we felt that a better understanding of the method in the absence of disorder was 
needed. 

By calculating exactly the two-point correlation functions for a random walk on a quad- 
rant of the lattice (the so-called corner rule 0), we will show that the procedure introduced 
by FG needs to be revised because of the uncontrolled approximations. This is also sup- 
ported by an exact analytical calculation on the one-dimensional analog. We will then 
proceed to show how one can obtain a consistent procedure based on the finite size scaling 



hypothesis 16 



The outline of this paper is as follows: In Sec.|TJ we first recall the derivation of the formal 
solution of the two-point correlation function for the random walk. This solution constitutes 
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the basis of the real-space renormalization group approach to this problem. The standard 
procedure is then reviewed and its fundamental problems are discussed. In Sec.^TT], the one- 
dimensional version of this problem is analyzed exactly in detail, and the inconsistencies of 
previous renormalization approaches are exposed. In SecJIV], we introduce our new method 
which is based entirely on the widely accepted finite size scaling approach. Then in Sec.[V] 
this method is applied to describe the phase diagram for a problem with different fugacities 
for the surface and bulk. The results are consistent with mean field theory where applicable. 
Finally Sec.[Vl| contains some conclusions and perspectives on the extension of this method 
to disordered systems. 



II. TWO-POINT CORRELATION FUNCTION FOR RANDOM WALKS 

It is well known that the unconstrained random walk on a lattice can be solved by 
using a generating function technique pioneered by Montroll [|1^,|1^,|1^]. Let P Xo ,x(.N) the 
probability for the walker to be at the position x G 7L d at the (discrete) time N, given that 
it started at the site Xq at the initial time 0. The master equation to be solved is then 

P Xo MN + l) = - E Px ,y(N) (2.1) 
y(x) 

where the notation y(x) means that the sum is restricted to the nearest neighbors of x. The 
number Cx ,x(N) = z N P Xo , x {N) of iV-step walks with end points Xq and x then satisfies 
the analogous equation, 

0^(^+1)= £ C X0! y(N) . (2.2) 
y(x) 

In order to solve the master equation, it proves convenient to introduce the generating 
function G XotX (k), 



G Xo ,x(k) = E k N C x<hX (N) 

E fc M (2.3) 



W.Xn^X 



where w is a walk having Xq and x as the end points, and |iu| is the number of steps 
associated with it. 

By multiplying ( |2.2| ) by k N+1 and summing over all N, one gets, taking into account the 
initial condition Ca; o ,a;(0) = 5 Xo ,x, 

(k) + 5x lh x ■ (2.4) 

y(x) 

It is easy to see that Gx ,x{k) is also the two-point correlation function of a scalar free-field 



theory and that Eq. ( |2.3| ) can be recovered from a von Neumann expansion (see e.g. [R 

Generally the procedure for a RG includes two basic steps: First one coarse grains 
microscopic details in real space or integrates over the fast modes in momentum space. 
This is followed by a rescaling of the space/momentum and of the model variables while 
conserving the partition function and recasting the Hamiltonian in the same functional form 
as before. 

In the case of the random walk, the requirement of the conservation of the partition 
function amounts to a mapping between the rescaled and original fugacities which can be 
formally expressed as: 

P({k'},w')=J2P({k},w) (2.5) 

w\w' 

where k' and w' are the fugacity and walk on the rescaled lattice, respectively, P stands 
for the partition function, and the sum is constrained to all w consistent with w'. In 
the case of self-avoiding walks, this procedure leads to a well known polynomial recursion 
relation between k' and k whose linearization around the fixed point leads to the value of 
the correlation length exponent v. On the other hand, once the self-avoidance is turned 
off, the polynomial recursion becomes an infinite series since there are an infinite number 
of walks even in the smallest possible cell. One is thus faced with the problem of finding 
either a way of summing over an infinite number of walks or a truncation procedure. Some 
time ago Family and Gould || devised a recipe along the latter line. Their idea was that if 
L = ba (a being the lattice constant) is the size of the system, then walks with number of 
steps iV larger than Nmax given by 



Nmax ~ [< Rn >}max ~ L 2 (2.6) 

will give a negligible contribution to the sum in Eq. ( |2.5| ). 

An L x L cell can be mapped into an L' x L' cell by the requirement Q2.5p as 

Xn(A;',^)=Xn(A;,^) (2.7) 

where we defined the quantity 

Xn(k,j)= Y. C n{^N)k N (2.8) 

and Cn(xo, N) is the number of A^-step walks in the cell starting from Xq and exiting in the 
direction h (n = 1, 2). According to the assumption of FG, x can be approximated as 

Xn(k,y)~ £ C A (x ,N)k N (2.9) 

L N=N MIN 

where Nmax is related to the system size by Eq. ( |2.6| ) and Nmin is the minimum number 
of steps needed to reach the closest exit of the L x L cell. 

As an example (cf. Fig.|l|), if we take V = a and L = 2a, then we have from Eqs. ( p.7|) 
and Q: 

A;' ps £ Cn(x ,N)k N = k 2 + 2k 3 + 5k 4 + Uk 5 . (2.10) 

N=N M in 

On the other hand the total number of A^-step walks having Xq,x as the end points can 
be calculated easily from Eq. ( p.4|) which gives the quantity ( |2.3|) exactly. Therefore both 
sides of Eq. (|2.7|) can be calculated exactly without any truncation procedure. In Fig.||] we 
compare the right-hand side of Eq. ( |2.8|) calculated exactly with the one calculated using the 
FG truncation procedure, which can clearly be recovered upon numerical Taylor expansion 
of the exact result up to the desired order. It is apparent that although in general the FG 
truncation procedure seems to reproduce rather well the trend of the fixed point k*, it fails 
to reproduce the singularity present in ( |2.8j ). This singularity moves closer and closer to the 
fixed point as the cell size increases. The physical origin of this singularity stems from the 
fact that, unlike other systems where criticality is reached only in the infinite volume limit, 
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the random walk has a criticality in any finite cell, by taking the limit N — > oo. This has a 
consequence, as seen in Table [I], that the value of the exponent v overshoots the exact value 
1/2 already at a very small cell size. This would also be the case with the FG truncation 
recipe if the size of the cell were pushed to a sufficiently large value (although the behavior of 
such an approximation scheme for very large cell size is not known and may be complicated, 
see 1). 

In the next section we will see how the same trend is found in the one-dimensional analog, 
where everything can be calculated analytically for any value of the system size L. 

III. THE ONE-DIMENSIONAL PROBLEM 

In order to gain firm understanding of the problems associated with the standard renor- 
malization approach, in this section we solve a one-dimensional version of the aforementioned 
corner rule renormalization method. Subsequently a semi-infinite one-dimensional model in 
the presence of an infinite, repulsive barrier at the origin will also be discussed. 

A. Transfer matrix solution for the one-dimensional corner rule 

Let us consider a one-dimensional lattice where x — 0, 1,2, ...,L and the sites x = — 1, 
L + 1 have an infinitely repulsive barrier (see Fig. |3|). (The lattice constant a is set equal to 
1 for simplicity.) The analog of Eq. (|2.4j ) for the correlation function for x > 1 is 

G , x (k) = k[G , x -i(k) + G , x+ i(k)] (3.1) 

along with the boundary conditions: 

G ,o(fc) = fcG ,i(fc) + l ( 3 - 2 ) 

and 

G 0tL (k) = kG ,L-i(k) , (3.3) 
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where we have assumed that all walks start from x = 0. 
We can put this equation in a transfer matrix form, 

* x (k) = T* x+1 (k) , x = l,2,..,L-l 



(3.4) 



where we defined 



1 \ 



V" 1 i J 



' G , x (k) ^ 

I Go,l-lW y 



(3.5) 



The eigenvalues of the matrix T are clearl 

l ± Vi - 4F 



A- 



2A- 



(3.6) 



which are real ifO<A;<l/2 and form a complex conjugate pair ifl/2</c<l (we consider 
< 1 in order to make the generating functions sensible). Note that A + • A_ = 1 in both 
cases. The right {U a } and left {U a } eigenvectors (a = ±), are given by: 



(l \ 



\ x± J 



/-a ^ 



U= 



A+-A. 



V 



1 



(3.7) 



where we normalized so that < U a \Up >= 8%. 

Then expanding ^i,(k) in terms of the right eigenvectors and using the boundary condi- 
tion ( p.3| ) we get the coefficients of U±, 



"A+-A. 



(3.8) 



Expressing ^i(k) by iteration of Eq. (|3.4j ) and using ( |3.8|) , we get 



Go,L , , l+1 \ L+l , 



A^ +i ) 



(3.9) 
(3.10) 



A+- A_ 

which determine Gq,l if we further impose the other boundary condition ( |3.2| ). Thus we 



find, 



A i - A_ 



A^ +1 - A^ +1 - k(\ L + - X L ) 



(3.11) 



for the system size L = 1, 2, .... 

Note that, although this solution is valid for arbitrarily large L, the boundary condition 
( p. 3D makes the system finite. This distinction will become more clear in the next subsection. 

Once Go^(k) has been obtained, the one-dimensional analog of Eq. (|2.8j), namely 
Xi(k, 1/L) = Go }X (k)\ x= L, gives the recursion relation corresponding to (|2.7|): 

Xt(V,l/L')=Xi(k,l/L) (3.12) 



The fixed point k*(L,L') is obtained by setting k = k' = k* in (|3.12|) . Using (|3.11|) and 



( |3. 12| ) for renormalization from a cell of size L + 1 to that of size L, we get an implicit 
solution for the fixed point k*(L + 1, L), 



\*L+2 _ \*L+2 _ I \*L+1 _ \*L+1> 

k* T = k*(L + 1.L) — r ^ T ^ { + T (3.13) 



where A^_ is \± evaluated at k = k* L . 

We expect k* L to approach 1/2 as L — > oo because the critical fugacity k c for the infinite, 
one-dimensional problem is exactly 1/2. Given this limit, we now look for the L dependence 
of k* L for large but finite L. For this purpose, we need to distinguish the following two cases 
due to the square root in the expression ( |3.6| ) for A: 

(a) Case < k* L < 1/2. 

If this were the case, the eigenvalues A?j_(/c) would be real. By introducing an auxiliary 



variable a by tanh a = y 1 — 4fc£ 2 and expressing ( p.!3|) in terms of hyperbolic functions, 
we get 

_ (1 + tanha) L+2 - (1 -tanha) i+2 - [(1 +tanha) i+1 - (1 - tanha) L+1 ] cosh a 
' ~ (1 + tanha) L+1 - (1 - tanha) L+1 - [(1 + tanha) L - (1 - tanha) L ] cosh a ^ ' ' 



We see that the right hand side of ( p.14 ) tends to 1 as a — > (i.e., k* L -> (§)-). This shows 



that k* L cannot be in this range at least for large L. 
(b) Case 1/2 < k* L < 1. 



S 



Here the eigenvalues A± form a complex conjugate pair. Again it is convenient to intro- 



duce a new variable 9 by tan6> = V 'Ak 2 — 1. Using 9, (|3.11|) can be expressed as 

r / 7 \ = sin0 

°' M ' sin[(L + 1)0] - ksin[L9] 

and the fixed point equation ( |3.13| ) as 

[(k + 1) -2cos#]sin# 



(3.15) 



tan(L^ 



(3.16) 

we can use a small #* 



L) cos 20- cos # + fc(l -cosfl) 

where #2 is the value of 9 at = Now in the limit k* — > (^ 
expansion. Since the right-hand side of (|3.16|) diverges as 9* — > 0, small values of #* are 
achieved only in the large L limit (as expected). Thus, the aforesaid limit is equivalent to 
the limit L — > oo. 

A 1/L expansion of ( |3.16| ) and use of the relation tan -1 x = n/2 — tan" 1 1/x yields, after 
some algebra, 



_ 7T 1 5n 1 257r 1 1 



This translates, in terms of the fixed point value, to 



(3.17) 



1 / 7r \2 1 



(3.18) 



The critical value fc c is obtained exactly from /c£ in the L —>■ oo limit as expected. 
However the exponent z/ does not have the correct limiting value (which is 1/2) as we show 
now. 

Linearizing the recursion relation (|3.12 ) around the fixed point, we get the eigenvalue 

dG , L (k)/dk . 



. dk' 

~ ~dk^ =k,( - L ' L,) 

dG 0>L {9)- 1 /d9 



dG 0iL '{k')/dk' k =w=k*{L,L>) 



dG^{9')-yd9'\JSll, 



L') 



(3.19) 



Using the expression (|3.15 ) and the large L expansions ( 3.1 7]) and ( 3.18 ) for renormalization 
from a cell of size L + 1 to one of size L, we get after some straightforward but lengthy 
algebra, 
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This serves as the denominator in ( |3.19| ), while the numerator must be calculated by sub- 
stituting L + 1 for L in Eq. (|3.15|) , differentiating as in (|3.19| ), and then substituting the 
expansions ( |3.f 7] ) and (|3.f 8|) . Thus we finally obtain 

AL = l + (2 + ^)^ + 0(-i) (3.21) 

and 

lnf(L + l)/Ll 1 

v = lim L \ . = — = 0.2238... 3.22 

l-oo In Al 2 + tt 2 /4 v ; 

This exactly calculated exponent v therefore badly overshoots the correct value v — 1/2 
in the same fashion as in the numerical evaluation of the previous section. Up to this 
point, however, it could still be a consequence of the finite size L imposed by the boundary 
condition ( |3.3|) . In the next subsection we will remove one of the boundaries and work with 
a semi-infinite one-dimensional system. 



B. Generating function solution for the semi- infinite system 

Let us again consider a one-dimensional lattice where x — 0, 1, 2, .... and the site x — — 1 
has an infinitely repulsive barrier (see Fig. |^), but there is no longer a barrier at the other 
end. For this calculation we will exploit a different method, namely the Laplace-Fourier 
method (see, e.g. second reference in ||). 

The master equation for the bulk is: 

P 0jX (N + 1) = i[P o ,,-i(A0 + P , x+ i(N)} (3.23) 

for x = 1, 2, ... and the boundary condition at the surface x = is: 

P 0t0 (N + l) = ~P 0tl {N) (3.24) 

We introduce the Fourier transform for the semi-infinite line by 
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f+ n da 

P(q, N) = Y, e-^Po^N) — P , X (N) = / ^-e^P(q, N) . (3.25) 
x=0 J-n 2tt 



Multiplying ( |3.23| ) by e %qx and summing from x — 1 to oo, and using the boundary condition 
O), we get 



P(q, N + 1) = P(q, N) cos q - ^e^P Q)0 (AA) . (3.26) 
We now Laplace transform in N by defining the generating function, 

oo oo 

G(q, A) = £ A 7V P(g, iV) = £ e-^G 0>x (A) (3.27) 

Ar=0 x=0 

where Gq iX (X) is the generating function for the probability Pq jX (N). Then from equation 
( P-26j ) and the initial condition Po tX (0) = 5o tX we get 

G ?, A = 2 ; . 3.28 

1 — A cos q 



This allows the determination of Go,o(A). Integrating ( |3.28 ) over the first Brillouin zone, we 
get 

G - (A) = iTwm (3 - 29) 

where we have defined the integrals 

ImW= rg*L- (3.30) 
J -it 2ti 1 — A cosg 

for x = 0, 1, 2, ... This integral can be easily computed as a contour integral 

4(A) = -/ dz- *L (3.31) 

7T JC:\z\=l A[Z — Z + ){Z — Z_) 

where 



*± = ^ • (3.32) 

It is worth mentioning that these poles are the same as the eigenvalues X± derived in 
the previous subsection, since the relation between the two fugacities is A = 2k. Again the 
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case < A < 1 and 1 < A < oo have to be distinguished. In the first case the poles z± lie 
on the real axis and z_ is always interior to the unit circle, while z + is always exterior. In 
the second case the poles z± are complex conjugate and lie on the unit circle symmetrically 
with respect to the real axis. The result for the integral is 

I Q (X)z x if < A < 1 

I*M = < \ ~ ( 3 - 33 ) 



1 *+- 



A z+— z- 



if 1 < A < oo 



where 



l/VT^P if < A < 1 
/o(A)H (3-34) 
if 1 < A < oo 



We now proceed to compute the G 0tX (X) in the two cases. From Eq. (|3.28|) and Q3.29|) , 
we find by inverse transforming in q, 

in the case < A < 1, and 

G ^X)=I X (X) (3.36) 

for the case 1 < A < oo. 

Consider first the case of 1 < A < oo. In this case, we introduce an angle 9 by 

z± = e ±ie , tan^ = VW^l . (3.37) 

Then, from (|3.33|) and ( p.36|) , we get 



Go.l(A)- 1 = . (3.38) 

cos v sin LB 

For large L, this expression changes sign rapidly as 9 (or A) varies, and thus unacceptable 
as the correlation function on physical grounds. 

Therefore, we now consider the case < A < 1. If we define, as before, an angle a from 

z± = e ±Q , tanha = Vl - A 2 , (3.39) 
12 



we find 

Gol(X)- 1 = e aL [tanha+ - 6 " ] . (3.40) 

2 cosh a 

We note that this is an exact result independent of any renormalization procedure. Now 
since the only singularity of G ,l(A) is at A = 1 (or a = 0) where 9 ° X L diverges, we must 
interpret A = 1 to be the critical point A c . Thus, it is also consistent with the mean field 
theory for semi-infinite systems where the transverse correlation at criticality behaves as the 
separation L goes to oo as L - ^ -1 ) for d dimensions. 

Now, applying the idea of cell renormalization discussed before to this result, we turn to 
the equation for the fixed point, 

GV(A*) = GV(A*) (3-41) 

where L' < L — 1 and L = 1, 2, .... This gives the fixed point at a* = or A* = 1 exactly 
independent of L or L', which is consistent with A c = 1. Turning to the critical exponent, 
we calculate the eigenvalue 



A 



(3.42) 



by linearizing the recursion relation around the fixed point A* = 1. Note that, since Go,l(A) 
is singular at A = 1, what we are attempting to do is to linearize around a singular fixed 
point. After some manipulation we find 

A(L,L') = ±±±, (3.43) 

which is clearly wrong, since it would yield v = 1 + 0(1 /L)\ 

Even though the renormalization approach based on Eq. ( |2.7| ) with ( |3.40| ) results in a 
nonsensical exponent value, the result for the two-point correlation function Go,l(A) itself 
is correct. Indeed it is easy to check that it reproduces the exact results for the critical 
exponents v and 71 (a surface exponent [TjJ) if they are calculated directly from it. More 
specifically one finds in the grand-canonical ensemble, 
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Xi(A) = £ Co,l(A) A ^T (1 - X)- 1 ' 2 (3.44) 

L=0 

which is the exact result for the surface exponent 71 = 1/2 (see next section) and 

eW= ^ L r°^ U ^ (1-A)- (3.45) 

which again gives the exact result v = 1/2. The latter result can also be read off directly 
from ( |3.40|) by noting that a ~ y/1 — A near A* = 1. 

So what is wrong with applying Eq. ( |2.7| ) to our Go,l(A)? In the present calculation, 
L is not the size of the system (which is infinite), but it refers only to a site of the semi- 
infinite, one-dimensional lattice. Therefore, the ideas of cell renormalization has no basis 
of application in this case. In fact, for general dimension d, the transverse correlation at 
criticality should behave as 

G o, i (A)U=A c ~ Z5 ^ r (3.46) 

for L large. Since T)± = 1 for random walks, such behaviour cannot be consistent with Eq. 
( p.7|) for any d other than 1. (We just saw that it does not work even for d = 1.) 
If we assume instead 

G ,l'(X') = (j) a G , L (X) (3-47) 

for an unknown exponent a, then this would lead to 

X*=l-0(j) , (3.48) 
O(t) , (3-49) 



a + 1 L' 

for renormalization from L + 1 to L. The choice of a = 1 then leads to the correct limiting 
values as L — >• 00. (This is also true for any ratio L/L' as shown in the Appendix.) The extra 
factor (L/L) a could be considered to correspond to the res caling step of the renormalization 
transformation. However, this choice of a is not in agreement with the intuitive guess of 
— d + 2 — 77 j_ for d — 1. Rather, it would correspond to —d + 2 — 77 where 7] = is the bulk 
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exponent. Thus it is not straightforward to repair this approach in a satisfactory way, even 
where exact correlation functions are available, not to mention that in most cases we do not 
have such luxury. 



IV. FINITE SIZE SCALING AND SURFACE CRITICAL BEHAVIOUR 

We learned from the previous calculations that, in the usual form, the cell renormalization 
procedure cannot be consistent. That is, it is not assured that, as the size of the cell increases, 
the results for the critical exponents improve and become exact in the limit of an infinite 
cell. If an improvement is attempted by using exact correlation functions, an inconsistency 
is again found stemming from the basic recursion relation ( |2.7|) itself. On the other hand 
one might still hope that, as the cell size grows, the approximate renormalization procedure 
of FG's truncation recipe could give better and better results. Unfortunately, however, this 
also appears not to be the case. Therefore, we need an alternative procedure which is both 
consistent in principle and workable in practice. 

The clue of where to start comes from the result ( |3.18| ) where we calculated how the 
fixed point was becoming exact in the L — > oo limit. This is indeed compatible with the 



finite size scaling hypothesis (see e.g |16| and references therein), 



fc c (L) = k c (oo) + A(^) 1 /" (4.1) 

where k c (oo) = 1/2 is the exact critical fugacity in the infinite lattice limit and v — 1/2 in 

this case. The idea, therefore, is that one can estimate k c {V) by looking at the divergences 
of the (bulk) susceptibility defined as 

1 oo 
l A l X €AN=0 

ill E E G Xo ,x(k) (4-2) 



where C(xq, N; A) is the number of iV-step walks starting from a point xq and entirely 
contained in the volume A and Gx ,x{k) is its generating function. The subscript B refers 
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to bulk in the sense that the endpoints Xq and x can be anywhere in volume A. Then, either 
by fixing the exact value of v = 1/2 one can calculate k c (oo) or by fixing the exact value of 
k c (oo) one can calculate the value of the exponent v. 

The results for the square lattice are shown in Fig.^ and Fig.|5|, and they are consistent 
with the expected values. Indeed a best fit for both cases gives A; c (oo) = 0.25 ± 0.01 
and 1/v = 1.94 ± 0.02 and in the case of v improves if we include more and more terms 
corresponding to larger cell sizes. 

The presence of the surface also changes the entropic critical exponents as it is well 
known . Indeed if the system is sufficiently large to make the distinction between surface 
and bulk sensible, one can decompose the total free energy Fi as: 

F L (Ak, h, h t , ~) i ^ >1 L d f B (Ak, h, i) + L^fsiAk, h, h u i) (4.3) 

where h and hi are the external fields associated with the bulk and the surface respectively 
and Ak = k — k c . By differentiating twice with respect to the proper external field and by 
using the finite size scaling ansatz, one gets the well known general results 

Xb{L) L> ^ tfl v 

UL) L> ^L^ (4.4) 



Xi,i(L) 



at the critical values Ak — 0, h — h' — 0. Here we have defined the local susceptibilities: 

l aJi l X £dAXeA 

Xi,i(k) = 4rr E E G Xo>x (k) (4.5) 

XoedAXedA 

where we mean by dA the boundaries of the volume A. 

Thus, our new method forgoes the usual corner rule renormalization per se, and instead, 
calculates various quantities associated with a finite cell and interprets them in terms of 
surface finite size scaling. From Fig.|6] one can see that this method very accurately repro- 
duces the values predicted by the mean field theory |TB|J2^], namely 7/V = 2, 71/^ = 1, 

71,1/^ = _1 - 
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V. PHASE DIAGRAM FOR SURFACE-BULK RANDOM WALKS 



As an application of the method just described, we now present the finite size scaling 
solution of the problem of the interplay between the bulk A (with fugacity k) and the 
surface dA (with fugacity hi) (see Fig.0), based on the exact calculation of the finite cell 
susceptibilities. 

Physically the possibility of changing the strength of the surface fugacity with respect to 
the strength of the bulk fugacity allows the surface to make up for the missing bonds. Clearly 
one expects that if k\ is sufficiently strong almost all walks lie on the surface, and then the 
critical point and the universality class should both change: when all the "interactions" in 
the bulk are zero, the walks are not allowed to stay in the bulk and we have the adsorbed 
phase. Since k± = 1/2 and k — 1/4 are the exact critical values corresponding to an infinite 
surface (a line in this case) and an infinite bulk, one then expects a qualitative phase diagram 
as shown in Fig.§. 

For this calculation we eliminated the corner and imposed a periodic boundary condition 
into the vertical direction 2, while the horizontal direction 1 is of size L and has free edges. 
The result is shown in Fig.|| It appears that there is a tricritical point (called special point) 
which is the intersection of three different lines (corresponding to three different second order 
phase transitions). Below the special point there is the ordinary transition, where the bulk 
and the surface undergo a transition at the same critical point. Above the special point 
there is a line of surface transitions, which take place if k\ is bigger then the special point 
ordinate, where the surface goes into an ordered state (where the susceptibility is singular) 
while the bulk is still disordered, as well as another line called the extraordinary line where 
the bulk also becomes ordered. 

Our estimate of the special point SP is at k = 0.25 ± 0.01, and ki = 0.35 ± 0.01, 
corresponding to a ratio ki/k = 1.40 ± 0.06. This is in reasonably good agreement with a 
simple mean-field argument which would predict the ratio of 4/3. The errors were estimated 
graphically. 
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The qualitative features contained in this phase diagram also appear in the case of the 
self-avoiding walks |2T[ and in percolation [^2[ . 

Quite similar features are found in the case of a defect (which is a d — 1-dimensional 
surface inserted into a d- dimensional bulk). The phase diagram obtained looks very similar 
to the one for the surface. The special point SP is found when ki/k = 1 again consistent 
with mean field arguments. 



VI. CONCLUSIONS 

In this paper we have given a detailed analysis of the difficulties associated with the usual 
cell renormalization approach to the random walk problem and presented an alternative 
method to calculate critical properties which does not suffer from similar difficulties. This 
new approach is shown to give results which improve as the size of the cell increases, unlike 
the previous approaches which are shown here not to have this essential feature, relying on an 
uncontrolled approximation. Our approach is based only on the finite size scaling hypothesis. 
Using this approach we have computed the full phase diagram of the effect of the surface 
fugacity having a different value from the bulk fugacity, and calculated the exponents 71 and 
71 1 as well v. All the results are consistent with simple mean field arguments as expected. 
Analogous results are obtained for the case of a defect. 

The real challenge now is to use this method for the situations where no simple mean field 
results are useful and the exponents are unknown. A timely example of such a case is the one 
where the substrate on which the random walks are constrained is disordered or otherwise 
self-similar. Although this problem is potentially important as it has essential features of 
transport through disordered media, an important materials problem, the results of previous 
investigations |||7|] are not reliable. This extension is not expected to be straightforward, 
however, since various numerical approaches [|]^,[l^,|14],[nj all agree on the result that the 
susceptibility singularity for an ideal chain or trapped ant (which can also be thought as the 
limit of the self-avoiding walk in absence of self-avoidance) is not a simple power law in the 
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presence of strongly correlated disorder, but rather an essential singularity. Further work 
will then be needed in order to implement this extension. 
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APPENDIX: EXACT DECIMATION FOR AN ARBITRARY SCALE FACTOR 



We will show in this Appendix that, unlike in the corner rule renormalization, an exact 
decimation procedure (which gives exact scaling relations) always gives the exact critical 
point and the exact critical exponent irrespective of the choice of the scaling ratio I = L/L' . 
Without loss of generality, we will do this for decimation from L to a — 1 (a being the lattice 
constant), with the rescaling ratio of L. 

It is not hard to convince oneself that a decimation of the L sites next to the origin gives 
rise to the recursion, 

and 

G , x , =x/L (k') = (1 - 2k 2 A L (k))G , x (k) (A2) 
where Al(&) is such that 

A « = rrki • < A3 > 

In the previous expressions ki refers to the fugacity k when the lattice constant is L. We 
can show by induction that k* L = \ (VL G N). Indeed if we assume k* L = k\ = | = k*, then 
we would have from Eq. ( |A1|) , 

A;2 +1 = fc*^(A;*) (A4) 

where 

1 _ 2k 2 Ar(k) 

D ^ = ^'W i-^il) !" ' (A5) 

Thus the claim is equivalent to showing that Di{k*) = 1, which is easy to derive using Eq. 



Next we will show, again by induction, that v = \ independent of L. From 
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k' L+1 = k' L D L (k) . 
Upon differentiation with respect to k, we have from Eq. 

dk' L+l _dk' L+ dD L (k) 

Since by hypothesis, 

-q^\* = ^MAx{L)[D L {k)\ k=k * + k L [ — ] 



after some manipulations we get, at the fixed point k* — 



2 ' 



D L (k*) = l, ^^ = 2(2L + 1) . 



Therefore, substituting in (\K*7\), we get 



dk 

which is what we wanted to prove. 



9k ' L+U ={L + 1) 2 = \max{L + 1) 
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FIGURES 

FIG. 1. Corner rule for the case of the smallest cell 2a x 2a. By symmetry the number of walks 
starting from the origin and exiting in the direction 1 or 2 are equal. 

FIG. 2. Comparison between the exact and approximate G\{k) = x^^V-^) m the cases 
L = 1,2,3. The three solid lines correspond to the exact evaluations, while the two dotted lines 
are the approximate results as discussed in text. The intersections of the the solid lines correspond 
to the exact fixed points k* = 0.3157, 0.2950, while the intersections of the dotted lines correspond 
to the Family-Gould fixed points k* FG = 0.3470, 0.3108 for a scaling from L = 2, 3 to V = 1, 2 
respectively. 

FIG. 3. One-dimensional case. In (a) the sites x = —I and x = L + 1 have an infinitely 
repulsive barrier, corresponding to the corner rule in d = 1. The second case (b) is a semi- infinite 
one-dimensional lattice with an infinitely repulsive barrier only at x = — 1. 

FIG. 4. Finite-size scaling result for the critical value fc c (oo), which corresponds to the true 
critical point for the square lattice. A best fit over all points gives k c (oo) = 0.25 ± 0.01, while the 
exact value is k c = 1/4. 

FIG. 5. Finite-size scaling result for the exponent v. A best fit over all points gives 
l/u = 1.94 ± 0.02, while the exact value is l/u = 2. 

FIG. 6. Evaluation of the bulk and surface susceptibilities in the finite-size scaling approach. 
The estimates for xb (0)> Xi (A) and xi,i (+) are obtained from the slopes of the lines shown in 
the log- log plots. The exact values according to the mean field calculation are n f/v = 2,~j\/v = \ 
and 7i,i/V = — 1. 

FIG. 7. Example of a surface-bulk problem where the surface fugacity k\ and the bulk fugacity 
k are shown. 
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FIG. 8. A sketch of the expected phase diagram. The lines shown correspond to the ordinary 
transition (O), the surface transition (S) and the extraordinary transition (E). Also shown is the 
tricritical point called the special transition (SP). 

FIG. 9. Computed phase diagram for the surface-bulk problem of Fig.^. The points shown have 
been calculated for different system sizes, L = 4 (0)> L = 10 (A), and L = 20 (O). The dotted 
and dashed lines correspond to slopes k\/k = 1,2 respectively. The special point is estimated to 
be at k = 0.25 ± 0.01, h = 0.35 ± 0.01. 
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TABLES 

TABLE I. Behaviour of v as functions of the cell size b = L/a where a is the lattice constant. 
The first two columns refer to the present work, while the second two refer to the results using the 
approximate recipe of Family and Gould ref ||. 



Scaling length b/b' k* v k* FG v FG 

2/1 0.3156 0.5438 0.3470 0.5853 

3/1 0.2950 0.5132 0.3108 0.5571 

3/2 0.2770 0.4441 0.2920 0.5129 

4/1 0.2825 0.4937 0.2926 0.5412 

4/3 0.2711 0.4485 0.2743 0.4868 

5/1 0.2745 0.4792 0.2838 0.5398 

5/4 0.2640 0.4351 0.2693 0.5148 
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